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1. Introduction 

Levy flights are stochastic processes characterized by the occurrence of extremely long 
jumps, so that their trajectories are not continuous anymore. The length of these 
jumps is distributed according to a Levy stable statistics with a power law tail and 
divergence of the second moment. This peculiar property strongly contradicts the 
ordinary Brownian motion, for which all the moments of the particle coordinate are 
finite. The presence of anomalous diffusion can be explained as a deviation of the real 
statistics of fluctuations from the Gaussian law, giving rise to the generalization of the 
Central Limit Theorem by Levy and Gnedenko [U [2] . The divergence of the Levy flights 
variance poses some problems as regards to the physical meaning of these processes. 
However, recently the relevance of the Levy motion appeared in many physical, natural 
and social complex systems. The Levy type statistics, in fact, is observed in various 
scientific areas, where scale invariance phenomena take place or can be suspected (see |3]- 
[6] and references therein). Levy flights are a special class of Markovian processes, 
therefore the Markovian analysis can be used to derive the generalized Kolmogorov 
equation directly from Langevin equation with Levy noise [7j. 

The problem of escape from a metastable state, first investigated by Kramers [8], 
is ubiquitous in almost all scientific areas (see, for example the reviews [9|, [10] and 
Ref. [H]). Since many stochastic processes do not obey the Central Limit Theorem, 
the corresponding Kramers escape behavior will differ. An interesting example is given 
by the a-stable noise-induced barrier crossing in long paleoclimatic time series |12j . 
Another new application is the escape from traps in optical or plasma systems [13]. The 
main tool to investigate the barrier crossing problem remains the first passage times 
technique. But for anomalous diffusion in the form of Levy flights this procedure meets 
with some difficulties. First of all, the fractional Fokker-Planck equation describing 
the Levy flights is integro-differential, and the conditions at absorbing and reflecting 
boundaries differ from those using for ordinary diffusion. Levy flights are characterized 
by the presence of long jumps, and, as a result, a particle can reach instantaneously a 
boundary from arbitrary position. One can mention some erroneous results obtained 
in [H] (see also the related correspondence [I5]), because author used the traditional 
conditions at two absorbing boundaries. There are a lot of numerical results regarding 
the different time characteristics of Levy flights, but obtaining the exact analytical 
results remains an open problem (see Ref. [3]). 

In this work, starting from the backward fractional Fokker-Planck equation we 
investigate the barrier crossing event in the presence of Levy noise, by focusing on the 
nonlinear relaxation time. The paper is organized as follows. In the following section we 
shortly review some recent results on barrier crossing problems with different approaches. 
In section 3.1, the generalized equations useful to calculate the nonlinear relaxation time 
(NRLT) are derived. In section 3.2 we give the exact expressions of NRLT for free Levy 
flights and for a cubic potential profile. Finally we draw the conclusions. 
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2. Barrier crossing 

The particle escape from a metastable state, and the first passage time probability 
density have been recently analyzed for Levy fiights in Refs. [3], [12], [16]- [26]. The main 
focus in these papers is to understand how the barrier crossing behavior, according to 
the Kramers law [8], is modified by the presence of the Levy noise. Here we discuss 
briefly some results on the barrier crossing events with Levy flights, recently obtained 
with different approaches. 

The main tools to investigate the barrier crossing problem for Levy flights 
are the flrst passage times, crossing times, arrival times and residence times. We 
should emphasize that the problem of mean flrst passage time (MFPT) meets with 
some difficulties because free Levy fiights represent a special class of discontinuous 
Markovian processes with infinite mean squared displacement. As already mentioned, 
the anomalous diffusion in the form of Levy flights, for a particle moving in a potential 
proflle U{x), is described by the fractional Fokker-Planck equation [6] for the probability 
density function W{x,t) 

dW d d°'W 

-^-IU'(.)W]^D^ (1) 

where the Riesz fractional derivative is deflned as 



1 /•+-/(a; + e) + /(x-e)-2/(a;) 



K{a) Jo e 



l+a 



d^. (2) 



and 



K (a) = — — , (3) 

^ ^ r(a + 1) sin (7ra/2) ' ^ ^ 

with T{z) the gamma function and < a < 2. Due to the integro-differential nature 
of the equation ([T]), we cannot apply the usual boundary conditions at the reflecting 
and absorbing barriers of the system investigated. The particle, in fact, can reach 
instantaneously the boundaries from any position. 

The numerical results for the flrst passage time of free Levy flights conflned in a 
finite interval were presented in Ref . [3] . There, the complexity of the first passage time 
statistics (mean first passage time and cumulative first passage time distribution) was 
elucidated together with a discussion of the proper setup of corresponding boundary 
conditions, that correctly yield the statistics of first passages for these non-Gaussian 
noises. In particular, it has been demonstrated by numerical studies that the use of 
the local boundary condition of vanishing probability fiux in the case of refiection, 
and vanishing probability in the case of absorbtion, valid for normal Brownian motion, 
no longer apply for Levy fiights. This in turn requires the use of nonlocal boundary 
conditions. Dybiec with co-authors found a nonmonotonic behavior of the MFPT as a 
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function of the Levy index a for two absorbing boundaries, with the maximum being 
assumed for a = 1, in contrast with a monotonic increase for reflecting and absorbing 
boundaries. 

According to the Kramers law, the probability distribution of the escape times from 
a potential well with the barrier of height Uo, has the exponential form 

p(t) = ^exp|-^| (4) 

with mean crossing time 

Tc = Cexp|^|, (5) 

where C is some positive prefactor and D is the noise intensity. The barrier crossing 
behavior of the classical Kramers problem was investigated, both numerically and 
analytically, in Refs. [3], [19]- [21], where the role of the stable nature of Levy flight 
processes on the barrier crossing event was analyzed. Authors considered Levy flights 
in a bistable potential U (x) by numerical solution of the Langevin equation associated 
to the fractional Fokker-Planck equation ([1]) 



X 



-U' {x) + e^") (t) , (6) 



where (t) is the symmetric Levy a-stable noise. It was shown that although the 
survival probability decays again exponentially as in Eq. (jlj), the mean escape time 
has a power-law dependence on the noise intensity D 

where the prefactor C and the exponent depend on the Levy index a. Using the 
Fourier transform of the Eq. ([T]) 

-ikU' (-i^]W -D\k\''W , (8) 



dt V 

the authors derived the mean escape rate for large values oil/ D in the case of Cauchy 
stable noise [a = 1) in the framework of the constant flux approximation across the 
barrier. The probability law and the mean value of the escape time from a potential 
well for all values of the Levy index a G (0, 2), in the limit of small Levy driving noise, 
were also determined in the paper [21] by purely probabilistic methods. The escape 
times have the same exponential distribution (jl]). The mean value depends on the noise 
intensity D, in accordance with Eq. with = 1, and the pre-factor C depends 
on a and the distance between the local extrema of the potential. 

The barrier crossing of a particle driven by symmetric Levy noise of index a and 
intensity D for three different generic types of potentials was numerically investigated 
in Ref. [21]. Specifically: (i) a bistable potential, (ii) a metastable potential, and (iii) 
a truncated harmonic potential, were considered. For the low noise intensity regime. 
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the result of Eq. ([7]) was recovered. As it was shown, the exponent /i(a) remains 
approximately constant, /x ~ 1 for < a < 2; at a = 2 the power-law form of 
changes into the exponential dependence ([5]). It exhibits a divergence-like behavior as a 
approaches 2. In this regime a monotonous increase of the escape time with increasing 
a (keeping the noise intensity D constant) was observed. For low noise intensities the 
escape time process corresponds to the barrier crossing by multiple Levy steps. For 
high noise intensities, the average escape time curves collapse into a single curve, for all 
values of a. At intermediate noise intensities, the escape time exhibits non-monotonic 
dependence on the index a, while still retains the exponential form of the escape time 
density. 

The first arrival time is an appropriate parameter to analyze the barrier crossing 
problem for Levy flights. If we insert in fractional Fokker-Planck equation ([T]) a (5-sink 
of strength q (t) in the origin we obtain the following equation for the non-normalized 
probability density function W (x, t) 

dW d d^W 

— = -JU'(.)W]^D^-,(t)S(.), (9) 

from which by integration over all space we may define the quantity 

gW = -^y Wix,t)dx, (10) 

which is the negative time derivative of the survival probability. According to 
definition (ITOl) . q [t) represents the probability density function of the first arrival time: 
once a random walker arrives at the sink it is annihilated. As it was shown in the 
paper [19] for free Levy flights {U (x) = 0), the first arrival time distribution has a 
heavy tail 

g(t)~t^/"-2 (11) 

with exponent depending on Levy index a (1 < a < 2) and differing from universal 
Sparre Andersen result [27t [28] for the probability density function of first passage time 
for arbitrary Markovian process 



In the Gaussian case (a = 2), the quantity (ITTI) is equivalent to the first passage 
time probability density (|T2l) . From a random walk perspective, this is due to the fact 
that individual steps are of the same increment, and the jump length statistics therefore 
ensures that the walker cannot hop across the sink in a long jump without actually 
hitting the sink and being absorbed. This behavior becomes drastically different for 
Levy jump length statistics: there, the particle can easily cross the sink in a long jump. 
Thus, before eventually being absorbed, it can pass by the sink location numerous 
times, and therefore the statistics of the first arrival will be different from that of the 
first passage. The result (fT2|) for Levy flights was also confirmed numerically in the 
paper [26] . 
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3. Nonlinear relaxation time with Levy flights 
3.1. General equations 

The nonlinear relaxation time technique is more suitable for analytical investigations 
of Levy flights temporal characteristics, because does not request a constraint on the 
boundary conditions. According to definition, the nonlinear relaxation time (NLRT) 
reads 

^., Io^[P{t,x,)-Pioo,xo)]dt 

P(0,xo)-P(oo,xo) ' ^'^^ 

where P (oo, xq) = lim P (t, xq) and 

/•L2 

P{t,xo)= / W {x,t\xo,0)dx (14) 

Jli 

represents the probability to find a particle in the interval (Li,L2) at the time t, if it 
starts from point x = xq. Let us use the same procedure as for calculating the first 
passage time probability density (see [29]). If the random process x (t) is Markovian, 
the probability density of transitions obeys the following backward Kolmogorov's 
equation [30] 

dW(x,t\xn,0) f , , , , 

^ 1 ' ^ =L+{xo)W{x,t xo,0) 15 

at 

with the initial condition 

W{x,0\xo,0) = 5{x-xo). (16) 

Here (xq) is the adjoint kinetic operator. After integration with respect to x from 
Li to L2 directly in Eq. (ITSil and taking into account Eq. iHM we arrive at 

^^^^ = L+(xo)P(t,xo). (17) 
The Eq. ( ITTll should be solved with the initial condition following from Eq. ( |T6l) 

P(0,Xo) = l(L„L,)(Xo), (18) 

where 1(Li,L2) (x) is indicator of the set (Li, L2). 

According to Eq. fITTI) I/+ (xq) P {00, Xq) = 0, and after integration of this equation 
with respect to t from to 00 we obtain (see Eq. f|T8l) ) 

L'^ (Xo) Q (Xo) = P (00, Xo) - l(Li,L2) (xo) , (19) 

where Q (xq) is the numerator of the expression (JT3il . i.e. 

poo 

Q(xo)= / [P(t,xo) -P(oo,xo)]c/t (20) 
Jo 
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Finally, in accordance with Eqs. fll3p and fl20p the nonlinear relaxation time can be 
calculated as 

1 - P(oo,xo) 

with xq G (Li, L2). Although Eqs. f|T9l) and fl2T|) are a useful tool to analyze the temporal 
characteristics of Levy flights in different potential profiles U{x), obtaining the exact 
analytical results for the generic a parameter, characterizing the anomalous diffusion, is 
one of the unsolved problems in this research area. Even for some particular potential 
profile, like the cubic one, to derive a general expression of the NLRT as a function 
of the Levy index a is a non trivial problem. In the next section we derive a general 
differential equation useful to calculate the NLRT for arbitrary Levy index and we find 
a closed expression for the case of Cauchy stable noise excitation (a = 1). 

3.2. Levy flights in a cubic potential 

The forward fractional Fokker-Planck equation for Levy flights in the potential profile 

U (x) reads 

ow (Ml 0) ^ a_ ^^^^^ ^, ^ i^ ^'"^;,;;;!"- °\ (22) 

where < a < 2. It is easily to find from Eq. (l22l) the expression for the adjoint kinetic 
operator 

i.(,„)_C'(,„)^ + Z,_£^. (23) 



Substituting Eq. fl23l) in Eq. f|T9|) we arrive at 

^d'^Qixo) , .dQ{xo) 

— — 7^ u [xo) — - — 

d\xo\ axo 



D-i^ - U' (xo) = P (00) - 1(^,^,) (xo) , (24) 



because the probability P (00, xo) does not depend on the initial position of the particles. 
The Fourier transform of Eq. (12^ gives 



Q(A;) -P'lfcrg(fc) 



^—ikL2 ^—ikLi 



'?(*:) = Tr- / Q{xo)e-""»dxo (26) 



where 



27T 

and we took into account that in accordance with Eqs. f|T^ and fl20|) Q (±C)o) = 
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Solving Eq. fl25l) and using the backward Fourier transform, we can calculate the 
nonlinear relaxation time as (see Eq. fl2T]) ) 

T (xo) = / Q (k) e^'^^dk, (27) 

1 - P(oo) J^^ 

where xq G (Li, L2). It is easily to check from Eq. (!26l) that Q {—k) = Q* {k). Dividing 
the integral in Eq. fl27j) on two parts for negative and positive variables k and using this 
relation we easily arrive at 

T {x,) = ^ e^'^^dk^ . (28) 

As a result, it is sufficient to solve Eq. fl25|) only for positive values of k 



d\ ( . d 



U" I—] - ikU' 



i- 



Q {k) - Dk"Q (k) 



dk J \ dk 

= P(oo)5(A:) + ^ ^ , (fc>0). (29) 

This is one of the main results of this paper. By solving Eq. fl29l) for a particular potential 
profile U{x), we are able to calculate the NLRT by using Eq. fl28|) . for a particle moving 
in that potential. However the general solution of this equation strictly depends on the 
functional form of the potential profile U{x), and not for all the potential profiles there 
is a solution of this equation. 

We now consider two cases: (a) a free anomalous diffusion, and (b) a cubic potential. 

(a) For free Levy fiights {U (x) = 0): P (00) = 0, and from Eq. ( !29l) we have 

^ p—ikLi p—ikL2 

After substitution of Eq. fl30|) in Eq. fl28|) and evaluation of the integral we find finally 
for the case < a < 1 

^fa)= tw'riW7r - pi) 

21^1 [a + 1) cos {Tta/ 1) 

As it is seen from Eq. ( !3T|) . the nonlinear relaxation time decreases monotonically with 
increasing the noise intensity D and has a maximum as a function of initial position xq 
in the middle point of the interval (Li,L2). For Levy index 1 < a < 2 the nonlinear 
relaxation time is infinite as for free Brownian motion (a = 2). 

(b) Levy fiights in a metastable cubic potential with a sink at x = +00 

U{x) = -^ + a^x. (32) 
Substituting this potential in Eq. (!29l) and taking into account that P (00) = we obtain 
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d 

dk 



2 dQik) 



dk 



+ {k^a^ - iDk^'^^) Q (k) 



-ikL2 g— ifcLi 



2n 



{k > 0). 



(33) 



To solve Eq. ( l33l) we introduce a new function R{k) = kQ{k). After substitution of 
this new function, Eq. (l33|l can be rearranged as 



d'^R{k) 



+ [a^ -iDk""-^) R{k) 



{k > 0). 



(34) 



dk^ ' ^" " " /"V-/ 27rfc 
It is quite difficult to find an analytical solution of this equation for arbitrary Levy index 
a. Thus, we limit our further considerations to the case a = 1. Substituting a = 1 in 
Eq. ([31]) and representing its right part in the form of integral, we arrive at 



-^+(a^-,D)fl(i) = -- 



L2 



-'^Uy (fc > 0). 



(35) 



This linear differential equation can be exactly solved, and, as a result, we find (for any 
A; > 0) the finite solution in the form 



Q{k) 



coe 



-f^k—i-yk 



+ 



L2 



iky^y 



2vr Ai y^ + iD- a? 



(fc>0). 



(36) 



where 



(5 = a 



7 = a 



1/4 



sm 



1/4 



COS 



1 ( D 

- arctan — 
2 

1 

- arctan 



D 



(37) 



Because of Q (0) < oo, the expression in curly brackets of Eq. (l36i) should be equal to 
zero, and we easily find the unknown constant cq 



Co = -: 



dy 



2vr Jl^ y'^ + iD-a^' 
Substitution of Eqs. ([36]) and §^ in Eq. ([28]) gives 



T (xq) = — Re 

TT I .10 



dk 



L2 ^-iky 



-jSk—i^k 



dy \ . 



(3^ 



(39) 



k Jli y^ + iD — a? 
After changing the order of integration and evaluation of the integral on k we arrive at 

dy 



1 f^^ ( D 

"<-'^J.. {2'" 



[A{xo,y)] +y2B{xo,y) 



(40) 
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where 



and 



{xo - y f 



B{xo, y) = arctan 



xo -7 



/3 



y2 



2 2 

y -a 



TT 



- sgn(xo - y) . 



(41) 



(42) 



In the following Fig. [T]we report the behavior of the nonlinear relaxation time T{xo), 
calculated by Eq. (140|) . as a function the initial position of the particle for different 
values of the noise intensity D, namely D = 0.07, 0.35, 1.0, 3.0, 5.0. 



X 




Figure 1. Nonlinear relaxation time, in arbitrary units (a. u.), as a function 
of the initial position xq, for five values of the noise intensity D, namely: 
D = 0.07, 0.35, 1.0, 3.0, 5.0. The values of the parameters are: a = 1, Li = —10 
and L2 = +10. 



The potential parameter a (see Eq. f l32p ) is a = 1, and the interval boundaries are 
Li = —10 and L2 = +10. The integration step used to calculate T{xo) from Eq. fHOj) is 
Ay = 10~^. For the initial position of the particle we focus on the range of values around 
the potential well, that is we consider xq G [—2, +1]. A monotonic decreasing behavior 
of the nonlinear relaxation time is shown. The NLRT decreases with initial positions 
moving from the left of the minimum {xq = — 1) towards the maximum [xq = +1) of the 
potential and with increasing noise intensity. An overlap of the different curves appears 
near the maximum of the potential. This behavior could be ascribed to the role of 
initial positions near the maximum. For initial positions that are close to the maximum 
of the potential (xq = 1) the height of the barrier to cross decreases considerably and 



Barrier crossing with Levy flights 



11 



the probability of the particle to fall back into the potential well increases. For the role 
of the initial conditions in barrier crossing, with Gaussian noise, see Refs. [TOl [TT]. 

In Fig. [2] we report the log-log plot of the behavior of the NLRT as a function of the 
noise intensity D, for three initial positions of the particle, namely: Xq = —2.0, —1.0, 0. 
As we can see the decreasing behavior of the NLRT with increasing noise intensity is 
recovered (see Ref. [21]). 



-0.9 




Figure 2. NLRT (a. u.) as a function of the noise intensity D, for three 
values of the initial position of the particle, namely: xq = —2.0,-1.0,0. The 
parameter values are the same of Fig. [TJ 



4. Conclusions 

In this paper we obtain the general differential equation useful to calculate the nonlinear 
relaxation time for a particle moving in a cubic potential and with an arbitrary Levy 
index a. For Cauchy noise {a = 1) we obtain the closed expression in quadrature of the 
NLRT as a function of the noise intensity, the initial position and the parameters of the 
potential. A monotonic behavior of the NLRT as a function of the initial position of 
the particle is obtained in this case. For free anomalous diffusion the NLRT decreases 
monotonically with the noise intensity as in the presence of the cubic potential. 
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